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Xf} • Abstract. The Riemann equation u t + uu x = 0, which describes a one-dimensional 

accelerationless perfect fluid, possesses solutions that typically develop shocks in a 
finite time. This equation is PT symmetric. A one-parameter PT-invariant complex 
deformation of this equation, Ut — iu(iu x ) e = (e real), is solved exactly using the 
method of characteristic strips, and it is shown that for real initial conditions, shocks 
cannot develop unless e is an odd integer. 
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1. Introduction 
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The concept of PT-symmetric quantum mechanics was introduced in 1998 pp. Since 
then, there have been many studies of the properties of such quantum theories 
[21 El HI El E]- In an effort to understand the underlying mathematical structure 
of PT-symmetric quantum mechanics, there have recently been a number of studies 
of classical differential equations that are symmetric under PT reflection. These 
studies include an examination of the complex solutions of the differential equations for 
the classical-mechanical systems associated with PT-symmetric quantum-mechanical 
systems [21 [71 El [9]. Additional studies have been made of various PT-symmetric 
classical-mechanical systems such as the PT-symmetric pendulum [10J, the PT- 
symmetric Lotka-Volterra equation [Tl] [12] . and the PT-symmetric Euler equations 
that describe free rigid-body rotation [12] . 

Many well-known nonlinear wave equations exhibit PT symmetry. For example, the 
Korteweg-de Vries and generalized Korteweg-de Vries equations [13], the Camassa-Holm 
equation |15j . the Sine-Gordon and Boussinesq equations [13] are all PT symmetric. 
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Once it is known that a differential equation is VT symmetric, it is possible to introduce 
a continuous parameter e to deform the equation into the complex domain in such a 
way as to preserve the VT symmetry. Such a deformation gives rise to large parametric 
classes of new equations whose properties exhibit an interesting dependence on e. For 
example, in recent papers a complex family of nonlinear wave equations obtained by 
deforming the Korteweg-de Vries equation was examined jTH [T7J . 

In this paper we study a simplified version of the complex class of Korteweg-de 
Vries equations considered in Ref. [T6j ; namely, the first-order nonlinear wave equations 
obtained by removing the dispersive third-order-derivative term. This class of equations 
is the PT-symmetric deformation of the Riemann equation. The virtue of studying this 
simpler class of equations is that they can be solved exactly and in closed form for all 
values of the deformation parameter e, as we show in this paper. This allows us to 
determine how the properties of the complex solution depend on e. 

The well known and heavily studied Riemann equation 

u t + uu x = (1) 

describes a one-dimensional incompressible accelerationless fluid. The velocity profile 
u(x,t) represents the velocity of the fluid at the point x and at time t. For the initial 
condition 

u(x,0) = f(x) (2) 

the exact solution to ([I]) can be written in implicit form as 

u = f(x — ut). (3) 

The Riemann equation (CQ) is symmetric under VT reflection, where V represents 
spatial reflection and T represents time reversal. (Note that under VT reflection u does 
not change sign because it is a velocity.) A simple way to deform the Riemann equation 
into the complex domain in such a way as to preserve its VT invariance is to introduce 
the real parameter e as follows: 

u t - iu(iu x ) e = 0. (4) 

Note that as in quantum mechanics, we assume that T is an antilinear operator that 
has the effect of reversing the sign of the complex number i. For the special value e = 1 
this equation reduces to the conventional Riemann equation (CQ). We will examine the 
solutions to (jl]) that arise from the initial condition in ([2]). 

Let us recall how to solve the standard Riemann equation in (CQ). The Riemann 
equation is a typical example of a quasilinear partial differential equation (one that is 
linear in the partial derivatives) that can be solved using the method of characteristics. 
To do so we parametrize the initial condition ([2]) in terms of the parameter r: 

t = 0, x = r, u = f(r). (5) 

Next, we construct the system of three coupled ordinary differential equations known 
as the characteristic differential equations: 

dt ^ dx du 

ds ' ds ' ds 
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In general, when the characteristics for a quasilinear equation are solved subject to 
the initial conditions imposed at s — 0, the resulting solution has the general form 
x = x(s,r), t = t(s,r), u = u(s,r). Carrying out this procedure for the case of the 
Riemann equation produces the solution 

t = s, x — f(r)s + r, u = f{r). (7) 

The last step is to eliminate the variables r and s in favor of the original variables x 
and t, which yields the solution u(x,t) of (JTJ) given implicitly in ([3]). 

It is well known that for most choices of the initial condition f(x) the solution 
u(x,t) eventually becomes a multiple-valued function of x. This happens because the 
characteristic curves in the x — t plane for the Riemann equation are the straight lines 
x = f(xo)t + xq, as we can see from (J7j). These lines are in general not all parallel 
because / is not a constant, and thus some of the characteristics may eventually cross. 
It is physically unacceptable for the velocity field of a fluid to be multiple-valued. To 
avoid multiple- valued solutions, it is conventional to require that when two characteristic 
lines cross, a shock develops. A shock is a jump discontinuity as a function of x in the 
fluid profile u(x, t). A solution u(x, t) containing a shock is not different iable at the 
position of the shock and it is called a weak solution [13J. 

The purpose of this paper is to show that deforming the Riemann equation into the 
complex domain has the effect of softening the shock singularities, as was conjectured in 
Ref. |14j . Specifically, we will show in the next section that shocks can only form when e 
takes odd-integer values. To show this we will solve the deformed Riemann equation (J4]) 
exactly using the method of characteristic strips, which is explained in the Appendix. 



2. Shock Formation Time for the PT-Deformed Riemann Equation 

To find out whether the PT-deformed Riemann equation (j3J) develops shocks, we 
examine the problem quantitatively: We calculate the precise amount of time required 
for a shock to form. To demonstrate the technique, we begin by calculating the shock- 
formation time for the standard Riemann equation in (CQ). 

A shock forms as soon as the coordinate transformation in ([7]) becomes singular. 
To determine when this happens we calculate the Jacobian of this transformation 



J = ^ = det 
d(s, r) 



l + sf(r). (8) 



1 

f(r) sf(r) + 1 

This transformation ceases to be invertible as soon as J vanishes, namely, at the minimal 

positive value of s = — l//'(r). Thus, the shock-formation time i s hock is given by 

positive min r „ . ... 
4hock = xeR [-V/'0*0]- (9) 

For e 1, (jlj) is nonlinear, and not quasilinear, and the method of characteristics 
does not apply. Instead, we shall solve (BJ using the method of characteristic strips, a 
generalization of the method of characteristics, which in some cases can be used to solve 
nonlinear first-order partial differential equations. 



Does the complex deformation of the Riemann equation exhibit shocks? 



4 



The method of characteristic strips is explained briefly in the Appendix. Using the 
notation in the Appendix, we rewrite (jlj) as F(x,y,u,p,q) = 0, where t = y, p — u x , 
q = u t , and 

F = q — iu{ip) e . (10) 

The procedure is to seek a coordinate transformation of the form x = x(s,r) and 
t = t(s, r) such that the five arguments of F satisfy the five characteristic-strip equations 
in fT42|) . namely, 

dx dt du dp dq 2 

Ts = ~ eq/p > dS = 1 > Ts = il ~ e)g ' Ts = pq/u > T s =q,u - 

Here we have substituted ( ITUl) into ( H2|) and also substituted F = in various places, 
after taking the appropriate derivatives. Note that for e = 1, the equations for x, t, and 
u agree with the characteristic differential equations ([6]) for the Riemann equation (TjJ). 

Remarkably, the coupled system of ordinary differential equations (TTTT) admits a 
closed explicit solution{jj] 

p = G(D - es)~ 1/e , q = E{D-es)- 1/e , 

w = £(D-es) (£ ~ 1)/e , t = s + H } x = —eEs/G + 1, (12) 

where D, E, G, H, and I are functions of the parameter r and are yet to be determined 
from the initial conditions (T5J). 

Next, as explained in the Appendix, we verify that ( Fl2l) . which solves the 
characteristic-strip equations (TTTT) . constitutes a solution of the partial differential 
equation (jlj). To do so we impose the condition that F(x,t,u,p,q) = on the five 
functions in (Tl2l) . In fact, this condition is nontrivial, and it translates into the condition 

[iGY = -i, (13) 

which determines G. [For irrational values of e the solution to this equation has infinitely 
many branches, and we choose the branch that tends to — 1 as e — > 1. This is the solution 
that corresponds to the solution of the Riemann equation (pQ) as e — > 1.] 

Next, we rewrite the initial conditions u(x, 0) = f(x) in parametric form. On the 
initial curve at s = we have 

x = r, t = 0, u = f(r). (14) 
Since u = u(x, t), we have u r = pf + q§ . Also f = f'(r) } f = 1, § = 0. Thus, 

u r = p(r) = f'(r). (15) 
Also, since F = 0, we deduce from ( flOl) that 

q(r) = iu(r)[ip(r)] e = if(r)[if'(r)] e . (16) 

|| It is not always possible to use the method of characteristic strips to solve a nonlinear partial 
differential equation because the technical assumption that is made in ([40)l may not hold. 
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We impose the initial conditions by setting s = in (TT21) . Then, using (TTBT) and (TTEj) . 
we obtain 

p(0) = f{r) = GD~ X I\ q(0) = if(r)[if'(r)] e = ED' 1 '*, 

u(0) = f(r) = ED^ 1 , t(0) = = H, x(0) = r = I. (17) 

From the first two equations in (JTTJ) we solve for D and E in terms of f(r), f'(r), 
and G. Substituting these expressions for D and E in the third equation in ffl7|) gives 
an equation for G that coincides with ( fl3i) . Hence, we obtain a self-consistent result. 
Substituting G into the expressions for D and E gives 

= —i[if'( r )Y e i E/G = -f(r)[if(r)Y-\ (18) 

From the last two equations in (JTTJ) we have 

I = r, H = 0, (19) 

and gathering our results from (Tl2l . (|T8|) . and ffT9l . we finally obtain 

£ = s , x = r + e/(r)[if'(r)] e ~ 1 s, u = f(r). (20) 

Surprisingly, even though (j2J) is nonlinear, we have solved it in closed form. [Note that 
for e = 1 this solution reduces to the solution of the Riemann equation (TJJ).] 

The next step in the application of the method of characteristic strips is to solve 
for r and s in terms of x and t, and to substitute the expression for r into u(r) to obtain 
the desired solution u(x,t). This can be done if and only if the Jacobian 



J=?H = det 



dt dt 

ds dr 

dx dx 

ds dr 



| = i+^{/w[^wr 1 } (2i) 



d(s, r) 
does not vanish. 

Clearly, for real f(r), that is, for real initial conditions, J never vanishes unless e 
is an odd integer. This establishes the principal result of this paper. When e is an odd 
integer, the Jacobian may vanish, and a shock will form at time t s hock given by 

positive min / — 1 \ 

tshock= rem {^iUmmF 1 })- (22) 

2.1. Case e = 2; 

For e = 2 (|1J) can be reduced to the Riemann equation for a pure-imaginary velocity 
field as follows: Consider (jlj) at e = 2 and multiply it by u. Define v = u 2 /2. 
The resulting equation is v t = —iv 2 x . Take the x derivative of both sides and define 
w = iv x = iuu x . Clearly, for real u(x, 0) = f(x), w(x, 0) = if(x)f'(x) is pure-imaginary 
and w satisfies the Riemann equation wt + 2ww x = 0. The Jacobian ( 12TI) in this case 
is J = 1 + 2is-^[f (x) f (x)]\ J can never vanish for real f(x), in accordance with our 
general result. 
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2.2. Case e = 3: 

According to our general result, shocks may form at e = 3. In this case ( [22]) reads 

1 



t 



shock 



positive mm 
r G 1R 



■ dr {f{r)\nr)f};- (23) 

As a concrete example, we take the initial velocity profile to be the Cauchy 
distribution f(x) = l/(x 2 + 1). For this case, (123|) becomes 



t 



shock 



positive min 
x e M 



[l + x 



2\6 



(24) 



_24x(l -4x 2 ) 

This function has two maxima and two minima, which are located at x — ±0.216 621 
and at x = ±0.769 392. The lower positive minimum occurs at x = 0.216 621, and at 
this value of x, it equals 0.311791. Thus, this initial condition evolves into a shock at 
t = 0.311791. 



Appendix: Concise Summary of the Method of Characteristic Strips 

In this Appendix we give a concise recipe for using the method of characteristic strips 
to solve nonlinear first-order partial differential equations. This Appendix is brief, and 
we remark that a complete discussion of the method of characteristic strips would be 
rather long and complicated (see, for example, Refs. [TBI [19]). 

We wish to solve the general first-order partial differential equation 

F(x,y,u,u x ,u y ) = (25) 

for the unknown function u(x, y) subject to the initial condition u(x, 0) = f(x). [Thus, 
u x (x, 0) = f'(x), while u y (x, 0) is determined by the equation F — at y — 0.] To this 
end we introduce two more unknowns, p and q and rewrite ([215]) as 

F(x,y,u,p,q) = 0, (26) 

subject to the constraints 

u = u(x,y), p = Ux =p(x,y), q = Uy = q(x,y). (27) 

In order to solve this equation in a manner similar to the method of characteristics, we 
seek a coordinate transformation of the form x = x(s,r), y = y(s,r), which will have 
the effect of reducing (125]) to a closed system of coupled ordinary differential equations 
in the variable s. The initial condition will be parametrized in terms of r. 

We can immediately write down eight differential equations associated with the 
partial differential equation f[2"5l) : 

P = u x , (28) 
q = u„, (29) 
Py = Qx, (30) 
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— = p — + q c }y (3i) 

ds ds ds ' 

dp dx dy ,g^N 
ds x ds ^ y ds ' 

dq dx dy 

■j- = + ( 33 ) 
as as as 

F x + F u p + F pPx + F q q x = 0, (34) 

F y + F u q + F pPy + F q q y = 0. (35) 

It is easy to check that f[3"4")) and fl33|) lead to = 0, which is consistent with (126"j) . 
Next, we use (1301) to rewrite (I34"I) as 

F a + F u p + F pPx + F qPy = 0, (36) 

and similarly, we rewrite (135]) as 

F y + F u q + F p q x + F q q y = 0. (37) 

Since our goal is to find expressions for the parametric derivatives 4^, and 

^ in terms of x, w, p, and g, we must try to eliminate all reference to the higher 
derivatives p x , p y , q x , and q y . We proceed as follows: Using (J32l) and ( J33l) . we eliminate 

p y and g x as 

ds x ds ) I ds ^ x \ds v ds J I ds 

Substituting these expressions into (1351) and ( 1361) we obtain 

+ F */| = „. (39) 

It is now clear that the only way to obtain a closed system of five coupled ordinary 
differential equations is to impose the additional requirement that 

F^ = F- (40) 
p ds q ds' [ ' 



Upon imposing (1401) . we see that ( 1351) and ( 1361) assume the particularly simple form 

F x + F u p + F q -£/-£ = 0, F y + F u q + F p ^J^ = 0. (41) 



We emphasize that the requirement (1401 is not necessarily valid, and after solving the 
partial differential equation (1251) it is necessary to verify the correctness of (1401) . 

Next, we make use of the functional degree of freedom to reparametrize; that is, 
to introduce a new parameter s' = g{s), which is an arbitrary function g(s) of the old 
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parameter s. This freedom allows us to choose ^ = F p , and from (1401) it follows that 



ds " P 

^| = F q . In this way we obtain the desired s derivatives as 



dx _ dy _ d« _ _ 

^ = -F x -F,p, %L = -F - F„g. (42) 
as as 

These are the characteristic-strip equations. These are a generalization of the three 

coupled characteristic ordinary differential equations for quasilinear partial differential 

equations to the five characteristic-strip equations for fully nonlinear partial differential 

equations. However, we emphasize that there is a new aspect; namely, that after solving 

the system PS]) , it is necessary to verify that the assumption made in (HUj) is valid. 

The problem of solving the nonlinear equation (1251) is now converted to solving the 

coupled system (HSl) . Once this is done and (|4TJ]) has been verified, it is then necessary to 

impose the initial conditions, and this will introduce the dependence on the parameter 

r. The final step is to eliminate r and s in favor of the original variables x and y, and 

thereby to obtain the solution u = u(x,y). 
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